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Abstract 

Many qualitative observations on the glass transition in classical fluids are well 
described by the mode coupling theory (MCT) but the extent of the non-ergodicity 
domain is often over-estimated by this theory. Making it more quantitative while 
preserving its microscopic nature remains thus a current challenge. We propose here 
a simple heuristic modification that achieves this for the long-time limit quantities 
by reducing the excess of static correlations that are presumably responsible for its 
inaccuracy. The location of the ideal glass transition predicted from this modified 
MCT compares very well with simulation for a wide range of interaction potentials in 
pure fluids and in mixtures. 

Being based on a time evolution equation deduced from a true hamiltonian, MCT 
is one of the few microscopic theories of the glass transition [1J. Developed initially for 
atomic fluids, its has more recently regained interest for describing the ideal glass transition 
in soft condensed matter [21 [3]. While it reproduces important qualitative observations 
in both physical situations, it suffers from quantitative limitations due to its approxi- 
mate treatment of the time evolution equation. Examples are the critical packing fraction 
rqMCT _ q_525 of the hard-sphere (HS) glass instead of if? ~ 0.58 from experiments on HS 
colloids [1], the critical temperature for the Lennard- Jones (LJ) fluid being in error 

by a factor two [5] and the whole T^ CT (p) curve for the square- well (SW) fluid [6] clearly 
misplaced in the temperature-density plane. Some ad-hoc recipes have been proposed to 
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correct this in specific situations such as the LJ mixture [7], but they are not completely 
satisfactory. Comparison with experiment or simulation is usually done by rescaling the 
state variables (see for example [6]), but at the expense of the microscopic nature of the 
theory 0[9]. To go beyond this, one should start from the basic approximations of MCT. 
Besides the projection of the dynamical variables on specific subsets, an important one 
is the factorization of four point contributions as products of pair terms [1] [2]. The static 
structure appears then in MCT at the level of the structure factors S a p{q) and the triplet 
direct correlation functions (dcfs) c a ^(q) (in Fourier space). The latter can usually be ne- 
glected (convolution approximation), as shown by previous estimates [TO] and confirmed 
by more recent ones [11] . As the liquid-glass transition is not very sensitive to the quality 
of the static structure, at least for hard-spheres [12J, it is tempting to view this overesti- 
mation by MCT of the non-ergodicity domain as reflecting too strong pair correlations, 
and attribute this to the factorization approximation. They should accordingly be re- 
duced (the possibility to predict the dynamics only from the static structure has recently 
been criticized |5]. However, the factored four point terms involving two time-dependent 
densities in q space, it is understandable that this approximation is more critical for the 
time dependence since different time scales might be improperly mixed). In this letter, 
we propose and test a pragmatic way of implementing this idea in the long-time limit: to 
reduce correlations, one should compute the static structure at an effective density that is 
lower (or a temperature higher) than the actual one (see [7] for a similar attempt). To see 
how this can be done in practice, recall that the densities p a = enter the MCT time 
evolution of the matrix S(q;t) with elements S a p(q;t) 



d_ 

dt 



S(q; t) + q 2 U(q)S- 1 (q)S(q; t) + J* dt'M(q; t - t')H-\q)-^S(q; t') = (1) 



through the irreducible collective memory function whose matrix elements are [1'6\ [T5] : 



55>{k;t) (2) 



M^(q,t) = 167r 3^ p ^(i/ 2 ) ]C J rfk ^;75(q' k )^;7'<5'(q> k )5'77'(llq- k ll;*)5'. 

In the vertex amplitude, 

V^ 5 (q, k) = h^M^C^k) + q.(q - k)5 M7 C^(||q - k||) + g\V 2 C^(k, q - k)] (3) 

the static dcfs are computed from standard methods of liquid state theory with the actual 
densities p a as input. The modification consists then in using effective densities p e J^ 
(specified below) in the vertex while retaining in the explicit dependence on density 
in the factor before the integral. As a result, the MCT equation is solved for the actual 
densities but with a static structure computed for p e J . Similarly, an effective temperature 
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T e f$ is used when temperature is a relevant variable. It is stressed that this is not a mere 
rescaling of the variables, a posteriori. Rather, this amounts to solving a modified MCT 
equation with an "effective " pair structure. To distinguish the results relative to the 
original MCT from those obtained from the modified MCT, we shall label the latter by a 
tilde (the label ex will refer to experiment, mostly computer ones here). 

To calibrate the correction, we start from the hard sphere potential for which there is 
a purely repulsive glass (caging mechanism). From the difference A 77 = 0.055 between the 
experimental critical glass packing fraction rjg X = 0.58 and the MCT one = 0.525 

with accurate static input [5j [11], we solve the MCT equation for the non-ergodicity 
parameter f(q) - at the packing fraction rj - by using the structure factor S(q;r] e ff) for 
an effective packing fraction rf^ = 77 — 0.055. The solution of this modified integral 
equation is the standard method of direct iterations [1] . The first non trivial result is that 
these iterations converge to a new critical value f^ CT = 0.584, virtually the exact one (ie 
7]g + A?]). Equipped with this calibration, we consider a binary hard-sphere mixture 
(size ratio 5 = D1/D2, packing fractions 771 and 772). The effective packing fractions are 
then 77? = rji — Ai](l — rjj/r}g X ); i ^ j = 1, 2. The slight dependence on the packing 
fraction of the other species is introduced so as to recover the one-component correction 
and prevent negative value. As with the one component case, we found that the modified 
MCT equations converge nearly to the corrected density, up to a size ratio 5 = 0.1. 
Quantitatively, the correction was tested on two HS mixtures considered in ref. |15| . For 
5 = 0.6 and x = rji/rj = 0.2, the total critical packing fractions are r]g 1CT = 0.528, 
~mct = 593i while ^ex _ 6Q6 from s i mmat i on . For 5 = 0.83 and x = 0.37, rj^ 101 ^ = 

0.524, fj^ = 0.589 and rf g x = 0.586. For both mixtures, the values predicted by the 
modified MCT are thus very close to simulation. This definitive improvement is confirmed 
by figure 1, which shows, as an example, the non ergodicity parameter fu(q) for 5 = 0.6. 
This excellent agreement with simulation - at the new critical density- is consistent with 
the relation rjg ~ 77^ — A77. It reflects the fact that the original MCT predicts the correct 
critical non-ergodicity parameter but a slightly inaccurate critical density, as found in [T2] . 

To go beyond the HS model, we first took the soft sphere potential 4> = e (f ) 12 f° r 
which density and temperature are combined in the coupling constant V = pa 3 (e/kBT) 1 ^ 4: . 
The one-component fluid of soft spheres has a critical coupling constant T e g x = 1.5 (see 
|16j). Using the Rogers- Young closure for the static structure as in ref. [T7], one gets 
Tg ICT = 1.33. Using as the only information the one gained from the HS potential, we 
compute the static structure in the modified MCT with an effective coupling constant 
f = T - 1.33(0.58/0.525 - 1). We then find a new critical coupling constant f g v ' ICT = 1.51 
again in very good agreement with simulation. We next took the mixture considered by 
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Barrat and Latz. One finds T g = 1.32 significantly lower than the simulation value 
pMCT _ i_4g ^3 j g re pi acea ; by = xfa^ + 2.x\X2o\ 2 + x\a\ 2 in the definition of T). 
With the modified MCT, the improvement is clear: T g ICT = 1.51. (using .575 instead of 
.58, we get T g ICT = 1.48 virtually the exact result). 

To consider a model with attractive contributions, we took the LJ fluid for which 
simulation data have been collected in [5] and more recently in [18] . Temperature plays 
now a role (the reduced temperature T* = ksT/e is used): a large error in the critical 
temperature Tg' ICT (p) can arise due to its very steep variation with p : a 10% variation 
of the density changes T g by a factor 2 [5[[7]. In this situation, slight details in the static 
structure might become relevant. The previous corrections for hard and soft spheres 
can be adapted by first noticing that the critical temperature T g {p) is well estimated 
from the critical packing fraction of hard spheres r\ g = ^pd? HS (T, p) where dns(T,p) is 
a suitably defined hard-sphere diameter (see eg. [5] and [2D]). In this spirit, we found 
(figure 2) that an even simpler formula ^pd? HS {T) = 0.58 fits very well the data of ref. 
[5] when the diameter is taken as the distance at which the reduced potential has a 
value ln(a) [19]: d H s(T) = <?[ 1+ j aT , ] 1/6 , with a = 0.8. T^ ICT {p) is also well fitted 
by a similar law with a = 0.5 for the effective diameter d e ff(T). Considering this as a 
seed, we computed the structure with an effective packing fraction rf^ * = rj — At/ with 
Ar ] = 0.58(l/d 3 HS (T)-l/d 3 eff (T)). 

The results is shown in figure 2. With the proposed modification, T g ICT (p) falls nearly 
exactly on the fitted law. As a check of the sensitivity to the choice of At/, the equation 
was solved by fixing it to the value for T g MCT (p) (ie the critical temperature predicted from 
the original MCT, for each density p). The results shown by triangles are close to the 
simulation of |18| which differ slightly from those of refs jo] (this difference is clearly visible 
on the enlarged scale used in figure 2). This might reflect different characterization of the 
glass (see e.g. [21]). The closeness of the triangles with the simulations of |18] is probably 
fortuitous, and there is a priori no reason to keep At/ fixed. With the definition of rf' ^ 
indicated above ( hence with a temperature dependent At/ ), we considered a standard 
LJ mixture at the state point studied in [7]: one finds T g ICT = 0.922 more than twice 
the simulation value T^ x = 0.435. With the modified MCT, we get T g MCT = 0.56, again 
definitively better than the original one (to speed the calculations, the static structure was 
computed from the closure of ref. [22]). 

Encouraged by these results, we finally considered the square well potential with short 
interaction range, a prototype for colloidal glasses for which an additional mechanism 
of dynamical arrest - the attractive glass - has been evidenced [3]. In this case, the 
MCT glass lines can be superimposed to simulation only with the help of a double linear 



4 



transformation of both temperature and density [6] . In figure 3, we show the result from the 
modified MCT equation using a temperature independent correction A77, a temperature 
dependent one Ary(T) with a linear dependence on 1/T* adjusted from a single simulation 
point besides the HS value, and lastly both Ary(T) and an effective temperature T*^ = 
T* + 0.18 ( to prevent cristallization, the simulation are actually for a slightly asymmetric 
mixture while we solved the MCT for the one-component fluid in order to compare with 
ref. |3]). Even in this more complex situation in which the mechanism for the arrest 
involves the formation of long-lived bonds, this way of reducing correlations leads to a 
clear improvement over the original MCT predictions. 

In conclusion, independently of the correctness of the view that attributes the quanti- 
tative insufficiency of MCT to the factorization approximation, there is little doubt that 
the method proposed here for tempering correlations definitely improves the accuracy of 
the present form of MCT, insofar as the long time results are concerned. Besides the 
operational value of this modification for accurately predicting the ideal glass transition 
for quite disparate potential, it is expected that these results that emphasize again the 
importance of many-body correlations and the presence of general mechanisms common 
to a wide class of models will stimulate studies from first principles. In the framework 
of microscopic theories, the extension of this idea to the time dependent quantities, by 
taking account different time scales, is for example conceivable. This should go in parallel 
with the development of more global approaches. Finally, this significant improvement of 
the MCT predictions that are relative to an idealized description of the glass transition 
should be useful in a step by step progress towards a better understanding of the complex 
phenomena that occur in real glass forming liquids. 
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Figure 1: Non-ergodicity parameter fn(q) in a binary hard-sphere mixture with D1/D2 = 
0.6 and x = 0.2 at the glass transition, Solid line: modified MCT (f]g ICT = 0.593); dashes: 
original MCT [12J CT = 0.5275); Symbols: simulation [E] = 0.606). 
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Figure 2: Glass transition line for the Lennard-Jones fluid. 
Filled squares and curve: original MCT and fitted curve with a = 0.5 in the effective 
diameter; Empty squares: modified MCT and curve adjusted to the simulation of [5] 
(filled circles) with a = 0.8; triangles: Modified MCT with fixed At/; empty circle: 

simulations of I IS I 
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Figure 3: predicted glass transition lines for the square well fluid. 
Dashes: original MCT; solid line without symbols: modified MCT with At/ = 0.055; 
filled circles: At/(T); squares : At/(T) and T e ^ . The theoretical results are for the 
one-component fluid with width 0.031c The dotted curve with empty circles shows the 

simulation data extracted from ref. [6]. 
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